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Abstract — System identification is a key enabling component 
for the implementation of new quantum technologies, including 
quantum control. In this paper we consider a large class of 
input-output systems, namely linear passive quantum systems, 
and study the following i dent inability question: if the system's 
Hamiltonian and coupling matrices are unknown, which of 
these dynamical parameters can be estimated by preparing 
appropriate input states and performing measurements on the 
output? 

The input-output mapping is explicitly given by the transfer 
function, which contains the maximum information about the 
system. We show that two minimal systems are indistinguishable 
(have the same transfer function) if and only if their Hamil- 
tonians and the coupling to the input fields are related by a 
unitary transformation. Furthermore, we provide a canonical 
parametrization of the equivalence classes of indistinguishable 
systems. For models depending on (possibly lower dimensional) 
unknown parameters, we give a practical identifiability condi- 
tion which is illustrated on several examples. In particular, we 
show that systems satisfying a certain Hamiltonian connectivity 
condition called "infecting", are completely identifiable. 

I. Introduction 

We are currently witnessing the beginning of a quantum 
engineering revolution [1], marking a shift from "classical 
devices" which are macroscopic systems described by deter- 
ministic or stochastic equations, to "quantum devices" which 
exploit fundamental properties of quantum mechanics, with 
applications ranging from computation to secure communi- 
cation and metrology [2], While control theory developed 
from the need for predictability in the behavior of "classical" 
dynamical systems, quantum filtering and quantum feedback 
control theory [3], [4], [5] deals with similar questions in the 
mathematical framework of quantum dynamical systems. 

System identification is an essential component of con- 
trol theory, which deals with the estimation of unknown 
dynamical parameters of input-output systems [6]. A similar 
task arises in the quantum set-up, and various aspects of the 
quantum system identification problem have been considered 
in the recent literature, cf. [7], [8], [9], [10], [11], [12], [13], 
[14], [15], [16], [17], [18] for a shortlist of recent results. 

In this paper, we focus on the class of passive linear 
quantum system which have been extensively studied in 
recent years [19], [20], [21], [22], [23], and can be used to 
describe or approximate a large number of physical systems 
[24]. The system consists of a number of quantum variables 
(e.g. the electromagnetic field inside an optical cavity), and 
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Fig. 1. Setup of system identification for linear quantum systems. The 
experimenter can prepare a time-dependent input state, and perform a 
continuous-time measurement on the output, from which the unknown 
system parameters 9 are estimated. The input-output relation is encoded 
in the transfer function H(s; 0). 



is coupled with the quantum stochastic input consisting of 
non-commuting noise processes (e.g. a laser impinging onto 
the cavity mirror). As a result of the quantum mechanical in- 
teraction between system and input, the latter is transformed 
into an output quantum signal which can be measured to 
produce a classical stochastic measurement process. In this 
context, we address the problem of identifying the linear 
system by appropriately choosing the state of its input and 
performing measurements on the output (see Figure [Q. 

In this "time-dependent input" scenario we characterize 
the equivalence class of systems with the same input-output 
relation (transfer function), find a canonical parametrization 
of the space of equivalence classes in terms of physical 
parameters, and investigate special parametric models in 
which the system is fully identifiable. These results are 
the first steps of an ongoing project aimed at designing 
statistically and computationally efficient input states, output 
measurements, and estimators, as well as characterising the 
class of equivalent systems in the "stationary input" scenario, 
which involves the analysis of the power spectrum. Although 
our investigations are guided and motivated by classical 
results going back to [25], [26], [27], the quantum setting 
imposes specific constraints and the goal is to identify 
quantum features related with e.g. unitarity, entanglement 
and uncertainty principles. 

The paper is structured as follows. In section II we in- 
troduce the set-up of linear quantum systems, and formulate 
the system identification problem. In section III we show that 
two minimal systems have the same transfer function if and 
only if their associated Hamiltonian and coupling matrices 
are related by a unitary transformation (cf. Theorem 12. IK 
amounting to a change of basis in the space of continuous 
variables modes. This yields a necessary and sufficient con- 
dition for the identifiability of a passive linear system, which 
is then applied to several examples with different Hamilto- 
nian parametrizations. In section IV we provide a concrete 
characterisation of the identifiable set of parameters with 



direct interpretation in terms of the coupling and Hamiltonian 
matrices. Finally, in section V we characterize a wide class 
of identifiable quantum linear networks, by employing the 
concept of "infection" introduced in [12], [14]. 

We use the following notations: for a matrix A = (ay), 
the symbols A\ A T represent its Hermitian conjugate and 
transpose of A, i.e., A^ = (a*J and A T = (aji), respec- 
tively. For a matrix of operators, A = (dy), we use the 
same notation, in which case a*j denotes the adjoint to dy. 

II. Linear quantum systems and identifiability 

In this section we briefly review the framework of linear 
classical and quantum dynamical systems and then formulate 
the quantum identifiability problem. 

A. Classical linear systems 

A classical linear system is described by the set of 
differential equations 



dx{t) 
dy(t) 



Ax(t)dt ■ 
Cx(t)dt 



Bu(t)dt, 
Du(t)dt, 



where x(t) € M n is the state of the system, u{f) £ K m is an 
input signal, and y(t) £ M. k is the output signal. The observer 
can control the input signal and observe the output, but does 
not have access to the internal state of the system. The input 
signal can be deterministic, in which case we deal with a set 
of ODEs or stochastic, in which case the equations should 
be interpreted as SDEs. Apart from the input and the initial 
state of the system, the dynamics is determined by the (real) 
matrices A, B, C, D. 

To find the relation between input and output it is conve- 
nient to work in the Laplace domain. The Laplace transform 
of x(t) is defined by 



C[x\{s) := / e~ st x(t)dt, 
Jo 



(1) 



where Re(s) < 0. Then, the following input-output relation 
holds: 

C[y](s)=E(s)C[u}(s), 



where 



S(s) = C(sI-A)~ 1 B + D 



(2) 



is the transfer function matrix. System identification deals 
with the problem of estimating the matrices A, B,C,D or 
certain parameters on which they depend, from the knowl- 
edge of the input and output processes. From ((TJ it is 
clear that the observer can at most determine the transfer 
function S(s) by preparing appropriate inputs and observing 
the output. 

The identifiability problem is closely related to the funda- 
mental system theory concepts of controllability and observ- 
ability. The system is controllable if for any states Xo,x\ 
and times to < t\ there exists a (piece-wise continuous) 
input u(t) such that the initial and final states are given by 
x(to) — xq and x(t\) — Xi, respectively. This turns out 
to be equivalent to the fact that the controllability matrix 
C = [B, AB, . . . , A n B] has full row rank. 



The system is observable if for any times to < t\, the 
initial state x(to) = xq can be determined from the history 
of the input and output on the time interval [io,ii]. This is 
in turn equivalent to the fact that the observability matrix 
O = [C T , (CA) T , (CA n ) T } T has full column rank. 

The importance of these concepts for identifiability stems 
from the fact that if the system is not controllable or 
observable then there exists a lower dimensional system 
with the same transfer function as the original one. The 
former can be obtained from the latter by separating it's 
coordinates via a canonical procedure called Kalman decom- 
position. Therefore, in system identification it is natural to 
restrict the attention to minimal systems, i.e. systems which 
are both controllable and observable. As noted above, by 
appropriately choosing the input signal u(t), the observer can 
effectively identify the transfer function S(s), while other 
independent parameters in the system matrices are not identi- 
fiable. The following theorem gives a precise characterisation 
of systems which are equivalent in the sense that they cannot 
be distinguished based on the input-output history. 

Theorem 2.1: Two minimal systems (A, B, C, D) and 
(A 1 , B', C , D') have the same transfer function H(s) if and 
only if they are related by a similarity transformation, i.e. 
there exists an invertible n x n matrix T such that 



D' = D. 



A' = TAT~ 1 , B' = TB, C' = CT~ X , 

B. Linear passive quantum system 

A general linear quantum system with n continuous vari- 
ables modes is described by the column vectors of creation 
operators a* := [d*,...,d*] T and annihilation operators 

a :- 



[di, . . . , d n \ satisfying the commutation relations 



Sijl. 



d. ( d* - d*d 4 = [di,d*] 
The system has a quadratic Hamiltonian of the form 
H = a f fia 

^11 • • • Qln 



(3) 



.,a_ 



a, 





di 




d„ 



(4) 



with Jlannxn complex Hermitian matrix, and is coupled 
to m bosonic quantum fields B(t) = [Bi(t), . . . , B m (t)] T 
whose algebraic properties are characterised by the commu- 
tation relations 

£;(«)] = min{s,t}S i:i i, 

or alternatively by 

[bi(t),b*(s)} = 8(t- 8)5^1. (5) 

where b(t) = [b\(t), . . . , b m (t)] T is the white noise operator 
formally defined as b(t) = dB(t)/dt. 

The coupling between system and field is described by the 
following set of operators (cy is a complex number): 



L = Ca = 



en 



Cln 





di 




d„ 



More precisely, the joint system-field evolution up to time t 
is given by the unitary operator U(t) satisfying the quantum 
stochastic differential equation [28] 

dU(t) = UBt(t)L - VdB(t) - ^Aadtj U(t), 

where 

A:=-iSl- -CftC. (6) 

This type of system is called "passive", because the operators 
do not involve the creation process such as a* a* in H and 
a* in L, representing a purely dissipative evolution. 
The Heisenberg evolution of the system operators is 

o(t) = U(t)*aU(t), 

which by differentiation gives the equation 

da(t) = Aa(t)dt - CUB(t). (7) 

Similarly, the output process 

B out (t) = U{t)*B{t)U(t) 

satisfies the differential equation 

dB ou \t) = Ca{t)dt + dB(t). (8) 

The Laplace transform £[x](s) of x(t) is defined as in 
(G3, for Re(s) < 0. As we will be assuming that the system 
is stable, the initial state of the system is irrelevant in the 
long time limit, and we can set its mean to zero (d(0)) = 0. 
In the Laplace domain the input-output relation is a simple 
multiplication 

C[b Ut ]( s ) =E(s)C[b](s), (9) 
where H(s) is the transfer function matrix: 

E(s) := I - C(sl - A)- l tf . (10) 

C. System identifiability 

Broadly speaking, by system identification we mean the 
estimation of the parameters 57 and C which completely 
characterize the linear quantum system described by (|7) 
and dH). This task can be analysed in various scenarios, 
depending on the experimenter's ability to prepare the state 
of the input and the initial state of the system, and the 
type of measurements used for extracting information about 
the dynamics. Here we restrict to a scenario which is most 
suitable for quantum control applications, and is a close 
extension of the classical set-up. More precisely, we assume 
that the experimenter can prepare input field in a coherent 
state whose mean vector has a desired time dependence 

(b(t))=p(t), 

and can perform standard (e.g. homodyne and heterodyne) 
measurements on the output, but does not have direct access 
to the system. 

By equation (0, the frequency domain output modes 
b (tbj) are related to their input correspondents b(iuj), 
through the linear transformations E(iui). Note that since 



b (iui) must satisfy canonical commutation relations simi- 
lar to (0, the matrix 3(iuj) must be unitary for all oj, as it can 
also be verified by direct computation [19]. It is clear that 
the experimenter can at most determine the transfer function, 
and this can be done by preparing appropriate inputs (e.g. 
sinusoids with a certain frequency oS) and observing the 
outputs. The transfer function matrix H(s) can be determined 
from the relation between input and output fields 

C[(b out )}( s )^E(s)C[(3]( S ). 

More generally, we can consider that the system matrices 
are not completely unknown but depend on an unknown 
parameter 9 € such that 

(n,c) = (n(0),c(0)), (ii) 

and S(s) = 3(s; ff) correspondingly. The task is then to 
estimate the unknown parameter 9 using the input and output 
relations (see Fig. [TJ, and we define the identifiability of the 
system as follows. 

Definition 2.1: The system parameter 9 is identifiable if 
H(s; 9) = H(s; 9') for all s implies 9 = 9'. 

III. Identifiability conditions 

In this section, we provide some basic necessary and 
sufficient conditions for the passive linear system and 
(0 to be identifiable, in the sense of Definition 12.11 

A. Observability, controllability and minimality 

The concepts of controllability and observability have a 
straightforward, though arguably non unique, extension to 
the quantum domain. The system defined by (0 and <[8j 
is controllable if for any times to < t\ and any means 
xq,x\ £ C™ there exists a coherent input which drives the 
initial coherent state |xo) into the final state |aci) over the 
time interval [to,ti]. This is equivalent to the fact that the 
controllability matrix 

C = -[(7 t ,AC7 t ,...,A n C7 t ] 

has full row rank. 

Similarly, the system is observable if for any times to < 
fi, an unknown initial coherent state of the system can be 
estimated from the history of (coherent) input and output 
over the time interval [£o,ii]- This is equivalent to the fact 
that the observability matrix 

= [C T ,(CA) T ,...,CA n ) T f 

has full row rank. As in the classical case, if the system 
is not controllable or observable then there exists a lower 
dimensional system with the same transfer function as the 
original one. Therefore, we focus on minimal quantum 
systems. The following lemma shows that in the case of 
passive quantum linear systems we need to check only one 
of the controllability and observability conditions to verify 
that the system is minimal. 

Lemma 3.1: For the quantum linear system ([7} and dHJ, 
the controllability and the observability conditions are equiv- 
alent. 



Proof: From the result of system theory, (A,C') 
controllability is equivalent to the following condition: 

yA = \y => yC^O. (12) 

This implies that the following condition holds: 

zAl = (iz z& ^ 0. (13) 

To prove this, suppose that there exists a vector z satisfying 
zA' = fiz and zC^ = 0. This implies that zCl = — i\iz 
and zC^C = 0, thus we have zA = z(~itt - C^C/2) = 
—fiz. But now z is contained in the kernel of C', this is 
contradiction to (lilt . The condition ( fT3] > is the iff condition 
for (j4J, Ct) controllability, thus equivalently (A, C) observ- 
ability. ■ 

We here provide an additional result for later use. 

Lemma 3.2: The quantum linear system (O and <[8J hav- 
ing the minimality property is stable. 

Proof: Because of the minimality, the system satisfies 
the condition (TOT l; hence z' is an eigenvector of A and fi* 
is the corresponding eigenvalue. Then the relation zA' = fiz 
readily leads to Rc(/i) = — ||Czt||/2||zt||, which is strictly 
negative due to zC' ^ 0. Therefore A is a Hurwitz matrix. 

■ 

B. The identifiability conditions 

As noted above, by appropriately choosing the input signal 
j3(t), the observer can effectively identify the transfer func- 
tion 5(s), while other independent parameters in the system 
matrices are not identifiable. The following theorem gives 
a precise characterization of systems which are equivalent 
in the sense that they cannot be distinguished based on the 
input-output relation. 

Theorem 3.1: Let (f2i,Ci) and (il 2 , C 2 ) be two linear 
systems as defined in and (O, and assume that both 
systems are minimal. Then they have the same transfer 
function if and only if there exists a unitary matrix T such 
that 

n 2 =TQiT\ C 2 =CiTK (14) 
Proof: It is well known in the system control theory 
that two minimal systems have equal transfer functions 

C x {sl - A!)- x Cl = C 2 (sl - A 2 )- 1 Cl 

(we here omit the trivial constant term /) if and only if there 
exists an invertible matrix T satisfying 

A 2 =TA 1 T-\ Cl=TC\, C 2 = CiT-\ (15) 

Note that T is not assumed to be unitary. Using the second 
and third conditions we have C%(T'T) — C\, which further 
gives [T^T, C\C\] = 0. Also, applying the second and third 
conditions to the first one, we have Q, 2 = TVliT^ 1 . Then, 
because f2^ is a Hermitian matrix, [T^T, Oi] = holds. 
Combining these two results we obtain [T*T, Ax] = 0. 
Therefore we have 

C x Ax = C7i(TtT)Ai = CMT^T), 

which means that the observability matrix O satisfies O — 
OT^T. Because of the assumption O is of full rank, thus T 



is unitary. Therefore the conditions ([TBI are reduced to (TT41 . 

■ 

For a parametrized model the identifiability condition is 
given by the following. 

Corollary 3.1: Let (£1(9), C(9)) be a minimal system 
with unknown parameter 6 E 0. Then 9 is identifiable if 
and only if 

fi(0') = Tfi(0)T + , C(9') = C(9)T^ 

implies 9 = 9' for any unitary matrix T. 

The above result can be interpreted as follows. In the 
classical set-up the matrix T corresponds to the coordinate 
transformation x 2 = Tx\ where Xi is the variable of the 
linear system. In our case, the transformation is a = Ta 
and the unitarity of T means that the canonical commutation 
relation (O is preserved. Note that if the system contained 
a classical component, T would not necessarily be unitary. 
Similarly, if the system is not passive, then one needs to 
consider both a and d* as coordinates, and corresponding 
doubled up transfer matrices [29]. 

In addition to the above corollary, we give another criterion 
for testing the identifiability. Note this result does not require 
the minimality of the system. 

Lemma 3.3: The parameter 9 is identifiable if and only if 

C(9)n(9) k C(9)' ! = C(9')VL(9') k C(9')\ Vfc (16) 

implies 9 = 9'. 

Proof: For simplicity let us denote C := C(9),C := 
C(9') and similarly for £1 and A. By expanding the equation 
H(s; 9) = H(s; 9') with respect to s and comparing their 
coefficients, we have CA k & = C'A' k C'^ for all k, and 
thus 

C( - iil - ic t C)' £ C t = C'( - iCl' - -C*C') k C*. 

This A;-th order polynomial is constructed by the linear 
combination of C[(C t C) p o fl^C* 1 with p + q = k, where o 
means the symmetrization, e.g., for k = 3 

(c^c*) 1 o n 2 = (dc)ti 1 + n^dn + ^(dcf. 

The relations (TfoT l can be proved by induction with respect 
to k, noting that for k = 0, 1 we have CC* = C'C't 
and respectively CflC^ = C'Vt'C'*, and the k step can be 
obtained from the expansion of CA k C^ = C'A' k C"^, by 
using all the previous steps. ■ 

C. Examples 

We will apply now the identifiable conditions to some 
quantum linear networks. For simplicity, we will assume that 
the coupling matrix C is known. 

Example 3.1: Consider a chain with three nodes, as de- 
picted in Fig. [2] (a). The system couples with a single 
input field (wavy lines) through the first node, such that 
L = \plka\ and the coupling matrix is 

C= [V2k, 0, 0]. (17) 



(a) (b) 
5XXX©-0-(D 




Fig. 2. Examples of quantum linear networks; (a) chain, (b) tree, and (c) 
ring. In each case, only the first node is coupled to the fields (wavy lines). 
The straight lines represent interactions between nodes. 



The system Hamiltonian has two terms (straight lines) de- 
scribing the interaction between neighboring nodes, and is 
given by 

H = 0i(a*a 2 + aia^) + 6*2(0203 + a 2 ag), 



hence 



fi(0) 



0i 

01 2 
0a 



where 6\ and 2 are unknown coupling constants which need 
to be estimated. 

To test the identifiability we apply Lemma 13.31 and com- 
pute 



CO(0)C t = 0, CO(0) 2 C t = 2k0 2 , 
Cn(6) 3 & = 0, OT(t9) 4 Ct = 2k0 2 (0 2 + e 2 2 ). 

Therefore Cfi(0) 2 Ct = CCl(9') 2 C^ implies 9\ = 0f, and 
Cn(9) i C = CQ.(6') A & implies 9\(9\ + 2 ) = 9' 2 (9' l 2 + 
9' 2 2 ). The two conditions together give 2 — 9' 2 , if 01 7^ 
0. One can verify that no additional relations arise from 
CQ(6) k C^ = Cfi(0') fe C* t with k > 4, and therefore the 
parameters 0i , 02 are identifiable up to a sign. 

For this system, we immediately see that the controlla- 
bility condition (thus the observability condition as well) is 
satisfied. Thus we can apply Corollary 13. II and find that the 
unitary matrices T\ = Diag(l, — 1, 1), T 2 = Diag(l, 1, —1), 
and T 3 = Diag(l, -1, -1) satisfy CT* = C and trans- 
form the system with parameter 9 = (61,62) into that 
with the system with parameter (— 6\, — 62), (— 0i, 62) and 
(01,-02), respectively. This gives an alternative proof of 
the fact that the equivalence class consists of the parameters 
{(±0i,±0 2 )}. 

A third route is to look directly at the transfer function 
given by 



E(s) 



and note that the poles give us enough information to 
determine both 2 and 9\. Note when 0i = (i.e., there is no 
connection between nodes 1 and 2), H(s) = (s — k)/(s + k), 
hence clearly 2 cannot be estimated. 



Example 3.2: The second example is a simple tree shown 
in Fig. 12(b). Node 1 is coupled to the input field and is con- 
nected to both nodes 2 and 3, thus the system Hamiltonian 
is of the form 



H = Aa 2 a 2 
hence we have 



i(al&2 + aia* 2 ) + 9 2 (ala 3 + aiaj), 



0(0) 



9 



(18) 



The C matrix is the same as in ( TFTI i. 

The additional Hamiltonian Aa;^ is necessary for the 
parameters 9\ and 02 to be identifiable, due to the following 
reason. If A = 0, we obtain C*O(0) fc C t = (9\ + 2 ) 2fc when 
k is an even number while it takes otherwise. Therefore, we 
can identify 2 + 9\ but not the individual components. This 
is due to the fact that when A = we cannot distinguish 
the nodes 2 and 3. To break this symmetry we introduce the 
term Aa 2 a 2 and we get 

coc* 1 = 0, cn 2 c f = 9\ + 9 2 2 

cn 3 ^ = A0 2 , cn 4 ^ = (9\ + 0^) 2 + 2 a 2 . 



This implies that all parameters are identifiable up to the 
signs of 0i and 02- 

Example 3.3: The last example is a ring-type network 
shown in Fig. [2] (c), where only the first node is accessible 
like the previous two cases. For this system the Hamiltonian 
is given by 

H = 9 1 (ala 2 + aia 2 ) + 9 2 (a* 1 a 3 , + 0103) 
+ 9s(a 2 ai + a2a\) + 04(0364 + 03^4), 

thus we have 

"0 01 2 

01 3 



n(9) 



02 
03 04 



<?4 




(19) 



while the C matrix is given by C = [V2k, 0, 0, 0]. 

Let us first check the minimality. A direct computation 
shows that the observability matrix O satisfies det(O) = 
4k 2 (0i0 3 + 2 04) 2 (0203 - 0104)- Hence, if the parameters 
satisfy det(O) ^ 0, the system is minimal. Now consider 
two set of parameters and 0' related through a certain 2x2 
unitary matrix U as follows 



U 



' 01 " 




' 0'i ' 


, u 


' 03 " 






. ^ _ 




. 8'2 _ 




. #4 _ 







(20) 



Then, by the unitary matrix T = Diag{l, U, 1} satisfies C 
CT and 










of 2 u T 







s 3 - ks 2 + (91 + 9l)s ~ k9\ 


TO(0)T t = 


U6 X 2 


O 


C/034 


= n(9') 


s 3 + ks 2 + (9 2 1 +61)s + k6 2 2 







9^U T 








Hence, from Theorem 13.11 the system with the parameter 
has the same transfer function as that of the system with 
0' defined by (l20l , Thus, this system is not (completely) 
identifiable. 



IV. Identification of the system matrices 

Let us consider the situation where the transfer function 
S(s) has actually been constructed from the input-output 
data. In the SISO (single input - single output) case, this 
will be a rational function of the form 

3(8) = 1 



-IS" 



CO 



a 



which is completely determined by the coefficients Cj, a 3 \ As 
it is well known in systems theory, the coefficient matrices 





1 







' " 








, Bo — 







1 







-ao -a\ 


—0>n-l 




1 



C = [Co, • • • ,Cn_i] (21) 

give a realization of 3(s), i.e., 5(s) = l + Co(s7— Ao) _1 i?o, 
and any other set of coefficient matrices having the same 
transfer function can be generated via the following similar- 
ity transformation: 

A = TA T~ 1 , B = TB , C = CoT" 1 . (22) 

However, the matrices (ETT i do not satisfy the constraints 
imposed on passive linear quantum systems. This means 
that, for general T, the transformation (l22t does not yield 
the set of coefficient matrices of quantum systems; e.g., the 
relation B = — C^ does not hold. Clearly, in this case, the 
system matrices (O, C) cannot be reconstructed. This issue 
is serious, because from the physics viewpoint we are often 
interested in those system matrices rather than the transfer 
function. Therefore, we need to find a specific T so that 
the coefficient matrices (l22l satisfy the constraints and that 
the system matrices can be reconstructed. In this section, we 
provide a concrete procedure to achieve this goal. 

A. Reconstruction of system matrices 

Let (Ao,Bo, Co) be constructed from the transfer function 
of the quantum system (|7]i and (O (note that now it is not 
limited to SISO). Then, for a certain matrix T, the matrices 
(1221 satisfy the constraints ©, which immediately yields 
A + A^ + C^C = 0, and B = -C*. These conditions are 
written in terms of (Ao, Bo, Co) as 

(T^T)Ao + Al(T^T) + C\Co = (23) 

and (T^T)Bo = — C\. Now the system is assumed to be 
minimal, thus Aq is Hurwitz from Lemma 13.21 This means 
that the Lyapunov equation d23l has a unique solution T^T > 
0. Accordingly, we have the diagonalization T^T = UoAUq, 
where A > is a diagonal matrix composed of eigenvalues 
of T^T and Uq the corresponding unitary matrix. Then, T 
is fully characterized by an arbitrary unitary matrix U as 

T = uVAUl (24) 

where y/X is a positive diagonal matrix satisfying (VA) 2 — 
A. This is the similarity transformation generating the equiv- 
alence class of quantum systems. In general, (l24l is not uni- 
tary, because (A 0l B , Co) are not necessarily the coefficient 



Unitary equivalence class 




(A ,B ,Co) 



(T A T , T B , CqTo : ) 



Equivalence class 
via all similarity transformation 

Fig. 3. Unitary equivalence class of system's coefficient matrices, which 
is generated from (Aq, Bq, Co). We denote To = x/AC/q. 



matrices of a quantum system. Also note that in a practical 
situation we need to check if d24T > satisfies the other condition 

(rtr)Bo = -cl 

Now we can reconstruct the system matrices. Again it 
follows from © that TAoT^ 1 = -itt- (CqT- 1 )^ {CqT' 1 ). 
Using d24l >. then we can represent ft as follows: 



il = UVLoU\ 



fin 



Similarly, from C = CqT -1 we have 

C = {CqUoVA^)^. 



(25) 



(26) 



These are the equivalence class of system matrices con- 
structed from (Ao, Bq,Cq), which are directly obtained from 
the transfer function, e.g. via (|2H in a SISO case. Note 
again that (l25l l. (l26b are the specific representations of (IT4l 
in Theorem 13. II see Fig. [3] 

B. Example 

Let us consider a two dimensional SISO system with only 
single node accessible and assume that the following transfer 
function has been experimentally obtained: 

Z(s) = l + ^^ , 

s z + axs + ao 

where a ,ai > and ci are real numbers. (As we will 
explain later, ci = —2a! is satisfied.) For this transfer 
function we take the typical realization ( |2T1 i; i.e., 



Ao 





-ao 



1 

-ax 



Bo 



Co = [0, ci]. 



Note that B = — Cq does not hold in general. With this 
choice, the Lyapunov equation ( 1231 has the following unique 
solution: 



2ax 



rt 





which is equal to A, and now Uq = I. Thus, the equivalence 
class of the system matrices are given by (EBI l and (l26l with 



On = 





-iJa~o~ 







ColVa - ^ [0,-, 



This means that the system is a two-node chain with internal 
coupling strength ^Joq and it couples to the outer field 
for measurement with strength \J2a\. In particular, when 
choosing U = [0, — 1 ; i, 0], we have 



ft = 







(1(1 




c = o], 



which have exactly the same forms as those in Example 3.1 
with specifically 6*2 = taken. Note that the condition 
{T^T)B = -Cq yields d = -2a x ; indeed this relation 
is satisfied for the two-node chain network, as easily seen 
by again setting 62 = in Example 3.1. 

C. Direct reconstruction of system matrices from the transfer 
function 

In Section IIV-AI we have shown that the system matrices 
can be reconstructed through typical realization methods 
employed in classical system theory. We here present another 
procedure that directly reconstructs the equivalent class of 
system matrices, using the specific constraints on the quan- 
tum systems. 

We begin with the simple SISO model where the coupling 
matrix is of the form C = {V0, 0, . . . , 0) with 9 > is an 
unknown parameter; that is, we assume that only a single 
node is accessible, as in the previous examples. However, 
we do not assume a specific structure on ft and write it as 



ft 



fl n E 



(27) 



E(s) = 1 



ifln 



E(s + ifl)- 1 E^ 



where ft is a Hermitian matrix with dimension n — 1. In this 
case, the transfer function (TlOb is given by 

6 
2 

Note again that we are assuming that H(s) is identifiable; 
that is, H(s) is a known complex rational function. The 
parameters are then reconstructed as follows. 

First, through a straightforward calculation we have 

*(1 - E(s)) 



1 + m n /s + l /2s + E(s 2 + isO)" 1 ^ 
which thus leads to 

9 = lim s(l-E(s)). 

s|— ¥00 

Next, with the use of the knowledge about 9 obtained above, 
we can identify On using the following equation: 



o n 



lim 



*gg(g) + 1) 
2(H( S ) - 1) 



Now, 9 and On have been obtained in addition to 5(s). 
This means that the function S(s) := E(sl + i£l)~ 1 E^ 
is known. We diagonalize as = VAV* with A = 
Diag{Ai,...,A„_i}. Then, E(s) = EV(sI - A)- 1 (EV)^ 
is of the form 

,_ , s + iXi 



where E[ is the i-th element of EV . This implies that A; can 
be detected by examining the function E(iui); that is, — iA, 
is the value on the imaginary axis such that E(iuj) diverges. 
Then, (assuming that ft has non-degenerate spectrum) we 
can further determine \El\ 2 from 

\E>f = ( S + l \ l )Z( S )\ s= _ ik . 

Lastly, let us express E[ as E[ = e^^E^ with real phases 
<j>i and define $ = Diag{0i, . . . , 4>n-i}- Then, d27l i can be 
written 

\E'\ 
A 



1 

Ve"^ 



fin 



1 











where \E'\ = [\E[ |, . . . , As shown above, the 

middle matrix can be completely identified from the transfer 
function 5(s). Therefore, all the eignevalues of fl can now 
be determined. Let us now summarize the result. 

Theorem 4.1: The equivalence class of systems having a 
given transfer function 5(s) is completely parametrised by 
the set of parameters (6, On, |_E^|, Aj) € R 2n , which are 
directly computed from 5(s) using the above procedure. In 
particular, the system parameter 9 and the eigenvalues of ft 
can be identified. 

To describe the general case, we assume that the m x n 
matrix C is of rank to; roughly speaking, this means that all 
the injected input fields couple with the system. Furthermore, 
we assume m < n; in this case, without loss of generality, 
C can be expressed as 

C=(C,0), 

with C a mx m full rank complex matrix. Correspondingly, 
we represent fl as in the same form ( |27] i, in which case O n 
is a ?7i x to Hermitian matrix. Then, as in the previous case 
we have 

= lim s(l - 3(»). 

s|— >oo 

This means that C can be represented in terms of a known 
strictly positive matrix Co an d an arbitrary unitary matrix U 
as C = CqU. Moreover, 



Un n U^ = lim l-iClil-Zis^Co + is! 



»|- 



2^0^*0' 



which means that On can be determined up to the unitary 
rotation by U. Now, we are given 

S(s) = UE(sI + m)- 1 (UE)\ 

which means that, from the same procedure as in the simple 
case, we can determine the eigenvalues of O and EiE* 
from S(s). Then, consequently, the eigenvalues of ft can 
be constructed. 

V. Network identification; the infection 

CONDITION 

In this section we treat a special class of models, which are 
identifiable thanks to the special structure of the Hamiltonian. 
This class was shown to be identifiable in different contexts 
in [12], [14]. 




(a) 



(b) 



(c) 



(d) 



(e) 



Fig. 4. Infection property. The colored node indicates that it is infected, 
and the arrow indicates that the infection occurs along that edge. Through 
the steps from (a) to (e), the network becomes infected. 



Let V be the set of vertices representing the modes of our 
continuous variables system. The interactions between the 
different modes are modeled by the set of edges £ over V: 

£ c V x V, 

so that two modes i and j interact if they are connected 
by an edge. More precisely, we assume that the matrix f2 
describing the system Hamiltonian is of the form 



(28) 



where u>ij are unknown real coefficients which make up the 
parameter 8 and e» = [0, • • • , 1, • • • , 0] T is the basis vector 
having zeros except the zth element. We further assume that 
the coupling between the system and the field is known and 
specified by the matrix C whose support is spanned by a set 
of basis vectors {e^ : i £ 1} for some set of vertices I, the 
restriction of C^C to this subspace being strictly positive. 

The crucial property we will require of I is that it 
is infecting for the graph (V,£), which can be defined 
sequentially by the following conditions (see Fig. 0): 

(i) At the beginning the vertices in I are infected; 

(ii) If an infected vertex has only one non-infected neigh- 
bour, the neighbour gets infected; 

(iii) After some interactions all nodes end up infected. 

Roughly speaking, this infection property means that the 
network is similar to a "chain", where the neighboring 
nodes are coupled to each other. Such a chain structure 
often appears in practical situations, and as shown in [30], 
it can be fully controlled by only accessing to its local 
subsystem. Also it is notable that in general a chain structure 



realizes fast spread of quantum information [31] and is 
thus suitable for e.g. distributing quantum entanglement. The 
result we present here is that such a useful network is always 
identifiable under a mild condition. 

Theorem 5.1: Let 0,(8) be given by Eq. d28l l. and assume 
that the support of C is spanned by {e^ : i G 1} with 
(1,V,£) having the infecting property. Then, if the system 
is minimal, £1(0) is identifiable. 

Proof: Recall that by Theorem 13.11 two parameters are 
in the same equivalence class if and only if there exists an 
n x n unitary matrix U such that 



Q(0 2 ) = USI{0{)U\ 



and 



c = cu. 



(29) 



(30) 



The latter condition implies [U, C^C] = and in particular 
U commutes with projection P onto the support of C^C so 
that 



U 











V 



(31) 



with V unitary on the orthogonal complement of the support 
of C. Let us write the Hamiltonian in the block form 
according to the partition J = 1 U I c : 



n(0) = 


( fin 


(0) 


n 12 (0) 




(0) 


o 22 (0) 


Then Eq. (f29b implies 


that 






On (02 


) = 


ill 





fi 12 (e 2 
n 22 (e 2 ) = 

The identity d32b means that 



i) = u>i,j(Q 2 ), i,jei. 



(32) 
(33) 
(34) 

(35) 



Furthermore, since 1 is infecting, there exists at least one 
vertex io € I which is connected to exactly one vertex jo € 
l c , so that rii 2 (0) can be written as 



n 12 (8) = 



T T 

w io ,3a e «o e j W *o Jo e *o e j 

T 

/ j / j -■uv-M-'-'i 



w%,j(d)(eie 



Eq. (|33T > then implies 

Ui o ,j o (0i)Vej o = u ioi j Q (6 2 )e jo , 

which means that ej is an eigenvector of V and U!i .j (9 2 ) = 
exp(i(f>o)u!i .j (di) for some phase <po. But since the coeffi- 
cients of 0(0) are assumed to be real, this implies that 



i G 1. 



(36) 



Additionally, since Ve JO = e_y , a decomposition of the form 
(|3T1 i holds with the identity block supported by the index set 
l'=lU{j }. 

The same argument can now be repeated for the set I', and 
by using the infecting property, all vertices will be eventually 
included in the growing set of indices, so that at the end we 



have f2(#i) = 0(02). Consequently, from Corollary 13.11 the 
system is identifiable. ■ 
From this result, we now readily see that the chain system 
studied in Example 3.1 in Section III C is identifiable, since 
clearly this system is infecting. On the other hand, the tree 
and the ring structural systems examined in the same section 
are not infecting, hence Theorem 15. 1 1 states nothing about the 
identifiability of these systems; in fact, as shown there, the 
tree system is identifiable, while the ring one is not. 

VI. Conclusion and future works 

In Theorem 13.11 we have shown that minimal passive 
linear systems with the same transfer function are related 
by unitary transformations. This is different from classical 
systems where the equivalence class is given by similarity 
transformations, the reason for that being the more rigid 
structure of the system equations ([7} and (|8) in the case 
of passive quantum systems. The characterisation of the 
equivalence class for general "active" linear systems remains 
an open problem. 

Another important issue which was addressed only briefly 
here is how to actually estimate the system parameters. More 
precisely, what kind of input state should be chosen and 
what measurement should be performed on the output. The 
performance of such a design of experiment can be measured 
by using statistical tools such as Fisher information and 
asymptotic normality [16]. 

Finally, the system identification problem can be con- 
sidered in a different setting, where the input fields are 
stationary (quantum noise) but have a non-trivial covariance 
matrix (squeezing). In this case the characterisation of the 
equivalence classes boils down to finding the systems with 
the same power spectral density, a problem which is well 
understood in the classical setting [25] but not yet addressed 
in the quantum domain. 

Some of these questions will be addressed in the future 
work [32]. 
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